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We study the importance of lattice refinement in achieving a successful inflationary era. We solve, 
in the continuum limit, the second order difference equation governing the quantum evolution in 
loop quantum cosmology, assuming both a fixed and a dynamically varying lattice in a suitable 
refinement model. We thus impose a constraint on the potential of a scalar field, so that the 
continuum approximation is not broken. Considering that such a scalar field could play the role of 
the inflaton, we obtain a second constraint on the inflationary potential so that there is consistency 
with the CMB data on large angular scales. For a m 2 </> 2 /2 inflationary model, we combine the two 
constraints on the inflaton potential to impose an upper limit on m, which is severely fine-tuned in 
■ the case of a fixed lattice. We thus conclude that lattice refinement is necessary to achieve a natural 

f***) ' inflationary model. 

Q \ PACS numbers: 04.60.Kz, 04.60.Pp, 98.80.Qc 
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I. INTRODUCTION 



Loop Quantum Gravity (LQG) provides a method of quantising gravity in a background independent, non- 
perturbative way. Quantum gravity is essential when curvature becomes large, as for example in the early stages of 
the evolution of the universe. Applying LQG in a cosmological context leads to Loop Quantum Cosmology (LQC) 
' which is a symmetry reduction of the infinite dimensional phase space of the full theory, allowing us to study certain 
O" 1 aspects of the theory analytically. In particular, it has been shown that classical big bang 0] and black hole singu- 
larities @ are removed in LQC, in a well defined manner. The discreteness of spatial geometry, a key element of the 
(3JT)j full theory, leads to successes in LQC which do not hold in the Wheeler-De Witt (WDW) quantum cosmology. 

In LQC, the quantum evolution is governed by a second order difference equation, rather than the second order 
| differential equation of the WDW theory. This is at the heart of LQC's success at removing singularities 0, Of 
equal importance is the continuum limit of these equations, which can lead to a behaviour which is qualitatively 
Q\ ' different from the classical cosmological evolution Q. As the universe becomes large and enters the semi-classical 
C**«- , regime, the WDW differential equation becomes a very good approximation to the difference equation of LQC. 

By deriving the effective equations of motion in the continuum limit, a semi-classical approach can be used to 
investigate possible phenomenological signatures of the theory. In particular, it was hoped that LQC would help 
\& • overcome the extreme fine tuning necessary to achieve successful inflation in general relativity Q. However, it has 
been shown 0] that semi-classical corrections are insufficient to alleviate this difficulty. There remains the possibility 
that effects due to the underlying discrete nature of LQC may come to inflation's aid. Here, we will explore this 
possibility and show that the discrete nature of the early universe naturally introduces constraints on a subsequent 
inflationary era. 

It has been suggested, by heuristic considerations, that for LQC to be considered physical the scale at which the 
discrete structure of space-time becomes significant must shrink as the universe expands [8| . This lattice refinement has 
■ been put on a more rigorous footing in Ref. @ . In particular, Ref. @ has analysed implications of lattice refinement 
for the semi-classical behaviour of the dynamical difference equations. Studying the effect of different choices for the 
lattice spacings on the solutions of the Hamiltonian constraint, the authors of Ref. || deduced the stability properties 
of the various refinement models. Combining the stability analysis with a confirmation of the appropriate classical 
limit, the suitable refinement models were constrained @. 

In what follows, in addition to discussing the effects of LQC on an inflationary era, we will show that lattice 
refinement is necessary to support a massive scalar field (in our context the inflaton) , and we will generalise to matter 
components with an arbitrary scale factor dependence. In particular, we will argue that lattice refinement renders a 
successful inflationary era more natural. 

In Section 2, we briefly discuss elements of LQC which we will later use. We are interested in isotropic models 
described in terms of one dynamical parameter. In Section 3, we quantise the theory and derive the Hamiltonian 
constraint. We then take its continuum limit, assuming that the regulating length (which is proportional to the length 
of the holonomy) is constant. We also discuss quantisation ambiguities. In Section 4, we investigate solutions of the 
full Hamiltonian constraint for large scales, which is indeed the WDW equation for a massive scalar field. In Section 
4, we discuss lattice refinement and we deduce the constraints imposed on the inflationary potential. In Section 5, we 
briefly discuss non-flat geometries. We round up our conclusions in Section 6. 
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II. ELEMENTS OF LOOP QUANTUM COSMOLOGY 

We deal with isotropic models described in terms of the scalar field a(t), thus there is only one kinematical degree of 
freedom. In connection variables, it is parametrised by the triad component p conjugate to the connection component 
c. Their relation to the metric variable a is [loj 

\p\ = a 2 and c = k + 7a , (1) 

where over-dot denotes derivatives w.r.t. cosmological time, 7 is the Barbero-Immirzi parameter, representing a 
quantum ambiguity parameter of the theory, and k depends on the intrinsic spatial curvature (k = 0, ±1). We will 
concentrate on flat universes, thus k = 0. The triad has an orientation, determined by the sign of p, and since the 
scale factor has no orientation, we used an absolute value in the above equation. The triad p incorporates information 
about the three-geometry through the scale factor a, while the connection c has information about the curvature, 
which in our case (flat and isotropic geometry) is contained in the extrinsic curvature proportional to d. The canonical 
variables c,p are related through 

fcp} = yV , (2) 
where re = 8irG and Vq is the volume of the elementary cell V adapted to the fiducial triad. Defining 

p = V^ 3 p and c = Vq^ 3 c , (3) 

with the triad component p determining the physical volume of the fiducial cell, and the connection component c 
determining the rate of change of the physical edge length of the fiducial cell, one obtains 

{c, P } = ^, (4) 

independent of the volume of the fiducial cell. 

To proceed in the quantum theory, one follows the procedure used in the full LQG theory. Thus, in LQC one takes 
e Vt)C / 2 , with j^o an arbitrary real number and p as the elementary classical variables, which have well-defined operator 
analogues (lO |. 

Using the Dirac bra-ket notation and setting e VoC / 2 = (c|/x), the eigenstates of p are the basis vectors \fj): 

*|M> = ^|M>. (5) 

Using the volume operator V = |p| 3 / 2 , representing the volume of the elementary cell V, with eigenvalues — 
(k 7 %|/6) 3 / 2 , one gets 0] 

%h(^) 3/ V (6) 

Notice that since we will be concerned with the large scale behaviour of the LQC equations, we will typically neglect 
the sign ambiguity that arises because the triads can have two different orientations. We thus avoid numerous factors 
of sgn(^) appearing in subsequent equations. 

To define the inverse volume operator, one has to trace over SU(2) valued holonomies. Since there is a freedom 
in choosing the irreducible representation to perform the trace, an ambiguity, labelled by J, arises. More precisely, 
the half- integer J stands for the spin of the representation. At this point, it is important to note that the same 
ambiguity arises in the gravitational part of the Hamiltonian constraint, since one has again to trace over SU(2) 
valued holonomies. Usually one quantises the gravitational part of the Hamiltonian constraint using the fundamental 
J = 1/2 representation, and the ambiguity is only investigated for the matter part. This issue is further discussed in 
Section 11 

Let us use the J =1/2 irreducible representation of SU(2). The inverse volume operator is diagonal in the \fj) basis 
and is given by [l(| 
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/3 _ yl/3 



|M) , (7) 



where /io is proportional to the length of the holonomy. At this point, we would like to stress that the regulating 
length fig is the crucial parameter in the quantum corrections. In Eq. ([7]) the eigenvalues are bounded and approach 
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zero near the classical singularity, in contrast to the classical case where the eigenvalues diverge at the singularity 
fi = 0. The eigenvalues reach their maximum at a characteristic scale equal to /ioi while at larger /i they approach 
the classical values and below /Uo they are suppressed compared to their classical analogues [|[. The inverse volume 
is cutoff for small volumes and, as it has been shown (llj . the generalisation to higher J will push the cutoff region 
to larger volumes. 



III. HAMILTONIAN CONSTRAINT 



The dynamics are determined completely by the Hamiltonian operator, the gravitational part of which, in the 
fundamental representation, is @], 



c 2i 

S K 2 fi,7 3 /iQ 



tr£^( 

ijk 



h i h j h i h- h k 



where 



hi = cos 



Hoc 



-1(7; sin 



(8) 



(9) 



denote the holonomies, with <j\ the Pauli spin matrices; the trace is taken over the SU(2) indices. We quantise this 
by setting 

= p\fJ,) and e«>e/2| M ) = e Mo ^|^) = + ■ (10) 
The action of the self-adjoint Hamiltonian constraint operator, fi g = (C g +C])/2, on the basis states, is 
3 



H g \fi) 



{ [R(n) + R(fi + 4/i )] | M + 4^ ) - 4R(/*)|/i) + [R(p) + R(p - 4/m)] |/i - 4/io)} , (11) 



where 



R(jjl) = (^n/ef 2 |m + Mo| 3/2 - |m - Mo| 3/2 



(12) 



We impose the constraint equation on the physical wave- functions which are those annihilated by the constraint 
operator, i.e. they have to fulfil the Wheeler-De Witt equation, since the classical expression must vanish. Such 
states can be expanded using the basis states as |$) = ^(0)1^), with summation over values of /i and where the 
dependence of the coefficients on <j> represents the matter degrees of freedom. Since the states |/x) are eigenstates of the 
triad operator, the coefficients <J M (^>) represent the state in the triad representation. Thus, quantising the Friedmann 
equation along the lines of the constraint in the full LQG theory one gets the following difference equation [l2[ 



V, 



V, 



p.— 3/io 



V, 



fi—b^o 



4k 2 7 j ^5 



(13) 



where the matter Hamiltonian 7^^ is assumed to act diagonally on the basis states with eigenvalue 7i^,(/i). Equation 
(I13p is indeed the quantum evolution (in internal time fi) equation. There is no continuous variable (the scale factor 
in classical cosmology), but a label \x with discrete steps. The wave- function ^(^>) depending on internal time fj, and 
matter fields <f> determines the dependence of matter fields on the evolution of the universe. 



Continuum limit 



We consider the continuum (i.e. /i 3> yuo) limit of the Hamiltonian constraint operator acting on the physical states, 
i.e. Wgl^). Expanding around small regulating length, fi , we get 



3 / n-/h\ 3/2 1/2 

2{— ^ 



j4_ 

24/r 



(a 3 -/? 3 ) + o(|) + hi. 



gher orders 



(14) 
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From Eqs. (fT3f and (fl4|) we obtain, in the small po limit, that 



2k 2 7 3 /i V 6 / 



2 2 / 3 \ 

^{^+4^(0) + «V-W</>)} - ^{^MoWO - 2*^(0) + *M-4„o(0)} + ( 



(15) 



The above second order difference equation distinguishes the components of the wave- functions in different lattices of 
spacing 4po- Assuming that \& does not vary much on scales of the order of 4/io (known as pre-classicality [Hj|), one 
can smoothly interpolate between the points on the discrete function ^(0) and approximate them by the continuous 
function ^[fj,,(f>). Under this assumption, the difference equation is be very well approximated by a differential 
equation for a continuous wave-function. As we will discuss later, it is this assumption that can break down at large 
scales, leading to deviations from the classical behaviour. 

Let us consider the parameter Mo entering in the regularisation of the Hamiltonian constraint to be a constant real 
number (the case of lattice refinement will be addressed in Section HV| . In other words, let us assume a fixed spatial 
lattice throughout the whole evolution of the universe. Expanding the wave-function ^(p ± 4/io, 4>) in Taylor series 
we obtain: 



*(^± 4^o,0) = *(A*,<£)± 



-(4W>) + 7T 



dp v ^ u/ ' 2 dp? 
Applying Eq. ([16]) in the difference equation Eq. (fl~5"|) we get 



(16/$ +0 Uo 



dp 3 



36 



k 2 j 3 H \ 6 / 



3/2 



+ ( M o) + o (pi) + ■■■ = -w (,j)V7^, 



(16) 



(17) 



where • • • denote higher order corrections in /i . Classically, the matter part of the Hamiltonian constraint for a 
massive scalar field is just, 



Ti<t> = n 



with momentum and potential V{4>). To quantise the above Hamiltonian constraint we use 



-ih- 



and 4>^(p, 4>) = 4>^{p, 4>) 



(18) 



(19) 



while the expressions for the volume and inverse volume operators, V and V 1 , are given by Eqs. ([6]) and {7}, 
respectively. Taking then the continuum limit, one obtains 



From Eqs. ((T7J) and l[2"0"]). the full Hamiltonian constraint, \ H e + H^j \H>) =0, reads 



or, equivalently, in terms of p with 



p = K"/h(i/6 



one obtains 



+ $ (VS»<» *)> ~ 3,-3/ 2 ~ + 0) + (mo) 



(20) 



(21) 



(22) 



(23) 



which is just a particular factor ordering of the Wheeler-DeWitt equation for a massive scalar field. 
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B. Ambiguities 



Thus far we have neglected any quantisation ambiguities, which can have significant consequences for the predictions 
of LQC 0, [H. As it was discussed in Section 2, there is a fundamental ambiguity arising from our choice of 
representation for the holonomies. It is possible to account for this ambiguity in the eigenvalues of the inverse volume 
operator [ll[, however constructing the gravitational part of the Hamiltonian constraint in non-fundamental (J ^ 1/2) 
representations, results in a higher order difference equation [8j. For J — 1 the difference equation has fundamental 
step size equal to 2fi , and it is of order 8 J = 8 [l2j. Being of higher order, spurious solutions may arise [12j suggesting 
that representations of higher order than the fundamental representation may not lead to the correct quantum theory. 
Whilst it may appear more natural to choose the same representation for both the matter and gravitational parts 
of the Hamiltonian constraint, there is no a priori reason this should be the case. Here we set J = 1/2 in the 
gravitational part of the constraint since it is the most tractable, whilst leaving it unspecified in the matter part, to 
remain as general as possible. 

The J ambiguity arises only in the inverse volume eigenvalues, which are given in general by [8[ 



K'jh^LJ (J + 1) (2 J + 1) ± 



V mV 2L/3 



where < L < 1 is a second quantisation ambiguity coming from the classical Poisson bracket 



{c, \p\ L } = g7/tLsgn(p)|p| 



L-l 



Taking the fi >> po limit, Eq. ([24]) implies 



i | QL-2)(3J 2 + 3J~lK 



hp 1 



O 



H- 3/ V> 



(24) 



(25) 



(26) 



There are also ambiguities arising from the classical matter Hamiltonian we choose to quantise 0], [T3|. Classically. 

1 



Ha, = K 



pi 

* a 3n 



2a 3(n+ 



is identical to Eq. (|18[) . however they will differ in the quantum regime. Equation (|26|) implies 

n(L-2)(3J 2 + 3J-l)A«o 



1 



(27) 



(28) 



Equations (f2"6")l and (|2"5)) imply that all quantum ambiguities are of order O (/^o//i 2 ) , thus they do not affect the results 
of Section HH3 



C. Wheeler-De Witt Solutions 



From Eq. (f2"3")l we see that for large scales the equation we want to solve is the standard Wheeler-De Witt equation, 



where 



Vp - gjf 2 ,0) + (V5*(P. 4)) + ap 3/2 V(^(p, 0) = , 



6 



At this point we approximate the dynamics of 4> to be such that 

V(<b)=p s -^V^ 



(29) 



(30) 



(31) 
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with V<j, a constant. For an inflaton field during slow-roll inflation, V(tf>) wconst., so we expect i5 « 3/2 in this 
era (which should be during the continuum epoch to be compatible with observations) . By separation of variables 
0) = T(p)*(0), we get 



+ (VP T W) + «P*^T(p) « 



(32) 



Equation (|32|) has solutions 



T(p) = Cip 1 / 4 



25 + 3 



26 + 3 



where J and Y are Bessel functions of the first and second kind, respectively, and C\ and Ci are integration constants. 
It is clear that for 5 > 1/2, the period of the solution decreases with increasing /i. Without lattice refinement (to be 
discussed in the next section) the discrete nature of the underlying lattice would eventually be unable to support the 
oscillations and the assumption of pre-classicality will break down, implying that the discrete nature of space time 
becomes significant on very large scales. 

We can see this rigorously by using the Bessel function expansions, 



Jp(x) 



2 ( Pit tt 

COS X 

,TX V 2 4 

2 . / /3tt tt 

sm x 

irx \ 2 4 



for large x, which give 
T(p) * 



26 + 3 



T7v2a^/V^ 



-(25+l)/8 



Ci cos x 



V3tt 



2(26 + 3) 



Ci sin x — 



%/3vr 



2(2^ + 3) 



(34) 



where x = 2^/2^^%(26 + 3)~V 2<5+3 )/ 4 . The zeros of T(p) are at 



tan x 



2(5 + 2) 4 



C?2 



(35) 



which gives the scale of the n -zero to be, 



(26 + 3)ir 



2V2a 



[C + n] 2S + 3 



with 



C = tan 



Ci 
C 2 



V3 



2(2(5 + 3) 



± 



1 

4 ' 



(36) 



(37) 



a constant. Note that we take the n zero of tan(x) to be at x = (2n ± l)ir/2. 

We are interested in the separation between two successive zeros in the limit of large n (large p) , which is given by, 

2tt 



Ap = lim (Pn+l - Pn) = 



2aJV« 



(38) 



For constant, it is clear that for a ^ (a = corresponds to the case where the wave-function is linear and hence 
has an infinite period), there will be a scale at which Ap is of the order of 2K"fhiio/3 (i.e. when Afi is of the order 
of 4/xq)) at which point the assumption that the wave- function is smooth on scales of the order of 4/xo breaks down 
and the semi-classical description is no longer valid. If we want the end of inflation to be describable using classical 
general relativity then it must end before this scale is reached. During slow-roll inflation the scalar field potential is 
approximately constant so we can set 6 « 3/2. 
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For Ap > 2K-fhfio/3, Eqs. (|22)l . (|30| and (|38| lead to the following constraint for the scale, a/, at the end of inflation: 

3-7T 

° f < 9 7W= ' (39) 

where we have also used that a = Ip] 1 ^ 2 = {n^h^^/Q) 1 ^ 2 . We can then calculate the maximum number of e-foldings, 
iV c ] = ln(af/ai), where a,i is the initial scale, which can occur within the classical era, given by 

oi > (2k 7 Vo/3) 1/2 . (40) 
Thus, from Eqs. (|39[) and (|4"0")) the number of c-foldings N c \ must satisfy the constraint 

N cl < In | , 3?r ] . (41) 

Setting [H ^to = 3V3/2 and 7 w 0.24, the above constraint, Eq. (|41| . on the number of e-foldings in units of H = 1 
reads 

AT cl « In [ 4Erl , (42) 



0'P1 . 



with lp\ = v8ttG. Thus, to have N c \ e-foldings of inflation with negligible quantum gravity effects, the potential must 
satisfy 

< ome-^n-f . (43) 

If we, conservatively, say that only half of inflation takes place in the classical era, the above constraint leads to 

F « lO- 28 ^- 4 . (44) 

Clearly, the above constraint becomes even tighter if we increase the number of e-foldings during the classical era. 
Considering lattice refinement we will see that this strong constraint on the value of the inflaton potential (during 
slow roll) becomes considerably less fine-tuned. 



IV. LATTICE REFINEMENT 



In the derivation of the Hamiltonian constraint we have assumed that the parameter fj,o is a constant real number. 
This parameter is related to the edge length used in holonomics. Since holonomics refer to embedded edges and loops, 
a constant fiQ corresponds to a fixed spatial lattice, so that only edges of a given coordinate length are used, and 
consequently the loop size remains the same as the total volume increases. In what follows we relax the assumption 
of a constant fig and we investigate the consequences, and thus the necessity, of a parameter jl being dependent on 
the phase space variable p, or equivalently, on /1. 

Allowing the length scale of the holonomies to vary dynamically, the form of the difference equation, Eq. (| 13[) . 
changes. Since the parameter /i determines the step-size of the difference equation, assuming the lattice size is 
growing, the step-size of the difference equation is not constant in the original triad variables. The exact form of the 
difference equation depends on the lattice refinement used. Here we will take the particular case of 

-* A(m) = /"oM~ 1/2 • (45) 

This type of lattice refinement is suggested by certain intuitive heuristic approaches such as noting that the minimum 
area used to regulate the holonomies should be a physical area [5j, or that the discrete step size of the difference 
equation should always be of the order of the Planck volume [8] . This choice also results in a significant simplification 
of the difference equation, compared to more general lattice refinement schemes. 

The basic operators are given by replacing fj, with jl. Upon quantisation we have Q 



(46) 
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which is no longer a simple shift operator since jl is a function of /i. If we change our basis to 

f d M 2 3/2 



(47) 



we find 

e - iA *|i/> =e-*"°*|i/) = |^ + Mo) ■ (48) 
Thus, the basis is a much more natural choice than The action of the volume operator on these basis states is 

3/2 

and the self-adjoint Hamiltonian constraint operator acts as 



(49) 



where 



9M 



1/2 



(50) 



C/ (i/) = \v + jtio| - |f - Mo| 



(51) 



Notice that since /i is dimensionless, so is za 

Following the same approach as previously, we expand |\&) = ^ u ((f>)\v} and the Hamiltonian constraint becomes, 



U{u + 4 Mo ) + U (i/)J „ + 4 M o (0) + 2|i/|U(i/) (i/) + - 



f + 4/i 
16m 3 , /6k 7 3 ^ 1/2 



V - 4:fl( 



U{v- 4 Mo ) + U (y) *„_4 M0 (0) 



(i>) *„ (0) 



9 V h 

Now we take the continuum (v ^> mo) limit of Eq. ([52"]) to get 



In terms of m Eq. ((53)) reads 



16 /6k 7 3 \ 1/2 



9 V ft 



= 



M 



-1/2 



d_ 

9/i 



M 



1/2 



(m 3/2 *(m, 



8 ^ 6k7 e 
3 



3\ 1/2 



H (m) * (j*,4>) + o (mo) 







(52) 



(53) 



(54) 



which is just a different factor ordering of the Wheeler-DeWitt equation, Eq. (|2"Tj) . Substituting in for Ti^ from 
Eq. ([20]), and for p from Eq. (J22J), we get 



-1/2 9 

V 1 Tr- 
ap 



p- 1/2 ^(p 3/2 ^(p,. 



- 48 P - 3/2 ^||^ + ^V(^^(p, t) + (mo 



= 



Thus, in the large scale limit the equation we have to solve reads 



-1/2 jL 
dp 



/3y(</>)// 2 *(p,0) = O 



where 



"-4 



(55) 



(56) 



(57) 



This equation is analogous to Eq. (|2"3"|) and thus we are following the same procedure as in Section IIII CI We first 
separate the wave-function <f>) into ^(p, <fr) = T(p)$(</>) and we then approximate the dynamics of the inflaton 
field by setting V{</>) = V^p 8 ' 3 ' 2 to get, 



-1/2 d 

P T 
ap 



1/2 ^(p 3/2 T(p) 



(3V 4>P s T(p) = 



(58) 
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The solutions of Eq. (158)) are 



t (p) = c lP ^ J* { li/MpVwA + c-zp-^ Y ( 

y ' 25+3 I 2(5 + 3 J \ 26 + 3 I 



Expanding Eq. ([59| . as we did in Section Hll C) we obtain 
where 



3tt 



2(25 + 3) 4 



T + C2 sin x 



3tt 



2(2(5 + 3) 4 



X = 4^^(2(5 + 3)- V 25+3)/4 , 
with /? defined in Eq. ([57)) . The scale of the n th zero is 



(26 + 3)tt 



[C- 



where 



C = tan 



-CA 1 



+ 1-1 



C 2 J tt 2(26 + 3) 2 4 
Calculating the separation between two successive zeros, we obtain 



Ap = lim (p n+1 - p„ ) = —== 



(l-25)/4 



For the continuum limit to be valid, the wave-function must vary slowly on scales of the order of p c = 
using p = K"fhp/6, we arrive at the constraint, 



3/2 

Ap > 4/i ( ^ ) p 



( Klh \ 



-1/2 



V 6 / 

Form Eqs. ([6"4")> and ([63)) we obtain the following constraint on V g 

27tt 2 



V A < 



(3-25)/2 



For slow-roll inflation, V(^>) must be approximately constant, thus 6 rs 3/2. As previously, we use pa — 
7 w 0.24, to find the constraint on the inflationary potential in units of h = 1, 

V((f>) <, 2.35 x 10^ 2 Zp 1 4 , 

where lp t = Gh. Clearly, this is a milder constraint than the one obtained for fixed lattices, Eq. (j43l) . 



(59) 

(60) 
(61) 

(62) 

(63) 

(64) 
4/i. Thus, 

(65) 

(66) 
3V3/2 and 

(67) 



V. OPEN/CLOSED (k = ±1) GEOMETRIES 

In the previous sections we looked at the k = case only. The reason for this is that there are several conceptual 
difficulties associated with the quantisation procedure in the other cases, particularly for k = — 1 (l4l . [l5j . These 
difficulties arise because vectors no longer commute and hence the four holonomies previously used to regulate the 
Hamiltonian no longer form a closed loop. There are methods of dealing with this (see Ref. Q and references therein) 
that lead to discrete evolution equations of the same form as Eqs. (|13[) and (|52p. but with an extra term involving 
the curvature. 

The correct quantisation should be able to reproduce the Wheeler-DeWitt equation including a curvature term in 
a suitable limit, 



d 2 S( P , 
dp 2 



+ kS([i,(j)) + &n s S(fi,(f)) = 



(68) 
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where S(fJ,,<ft) — ^/^^(fi, <fi). This no longer has general solutions, however during inflation when j « 1, analytic 
solutions do exist, 

T(m) = CiAiry A (- (k + an) a" 2 / 3 ) + C 2 Airy B (- (fc + a/x) a~ 2 / 3 ) , (69) 
where T(/z) is the \x dependent part of S(n,4>) (i.e. S(n,4>) = T(/i)$(</>)). Using the Airy functions expansions, 

lim Airy A {-\x\) = J~ U4 sin ( \x 2 l z -J- - 

lim Airy B (—\x\) - 



'TTX 1 / 4 V 3 4 

1 ^2 2/3 77 



we obtain an explicit form for T(/i), which has zeros at the scales 



Mn=a- 1 / 3 ^J (C + nr 73 -^, (70) 
where now 

C = tan" 1 (-C2/C1) - 3/4 , (71) 

is again constant. Clearly then, the presence of a curvature will not affect the large scale behaviour of the scale 
between subsequent zeros and so we expect the previous section to hold for any quantisation schemes that have a 
continuum limit of the form of Eq. (|68[) . As we have seen, changing the factor ordering of the Wheeler-DeWitt 
equation often means that analytical solutions are not available, in particular adding a k^/p Y(p) term to Eq. (|56[) 
results in an equation that cannot be solved exactly. However, the success of this factor ordering (which arises from 
a non self-adjoint Hamiltonian constraint) suggests that the results for other cases may be similar. 

VI. DISCUSSION/CONCLUSIONS 

We have studied phenomenological aspects of LQC. In particular, we have investigated the effect of lattice refine- 
ment, which as we will now argue, makes the conditions for successful inflation more natural. 

Let us consider that the scalar field <f> plays the role of the inflaton. We will then combine the condition imposed 
earlier on the potential V(<f>), with the constraint so that the perturbations induced by <f> are consistent with the 
CMB measurements on large angular scales. As a concrete example we will use a simple inflationary model, namely 
V{</>) = m 2 <f> 2 /2. 

The contribution of the inflaton field <f> with potential V(4>) to the fractional over-density in Fourier space, at the 
epoch of horizon crossing, k = aH , is given by 



2/3 



75w 2 M^ [V'((b)Y 



(72) 



iH 



where V'(4>) = dV/d<f> and the Planck mass is M Pi = (SttG)" 1 / 2 ~ 2.4 x 10 18 GeV. From the COBE-DMR measure- 
ments, 5 H ~ 1-91 x 10~ 5 PI, which implies from Eq. (J72J) that 

« 5.2 x lO^Af 3 , . (73) 

To ensure the continuum approximation is not broken at large scales, we found in previous sections (Sections III 
and IV) that the potential must be V((f>) ^ x, where x depends on whether we have static or dynamically varying 
holonomy length scales. In the case of fixed lattice and assuming 7V c i e-foldings of inflation within the classical era, we 
found x fixod lattico - 7 x lO^e- 2 ^- 4 . Whilst, considering lattice refinement we got a; lattico rcfincmcnt (^) < lO" 2 /" 4 . 

We now consider the simple inflationary case V((f>) — m 2 (j) 2 /2, favoured by the latest WMAP measurements. In 
this case, Eq. (|73| implies 

m(j) 2 « 1.5 x 10- 3 Af^ , (74) 
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and the constraint on the potential so that the continuum approximation is valid on large scales reads 

m 2 cj) 2 < 2x . (75) 

Equations (73]) and (75]) imply 

=> to ^ lO^/Mlj . (76) 

For the fixed and varying lattices, the above constraint reads 

to < 70{e- 2Nd ) Mpi (77) 
and to < 10 A/pi , 

respectively. Clearly, for any significant proportion of inflation taking place in the classical era, Eq. (|77h ) is a very 
strong, and therefore fine-tuned, constraint on the mass of the inflaton field. The condition on the inflaton mass 
becomes however natural, Eq. (|77b ). once lattice refinement is taken into account. 

Equations (TTj) can be compared to restrictions on the inflaton mass coming from the WMAP data [l9j, for the 
same inflationary model as the one we consider here, namely: 

/ <rn \ 

- 5.24 < log < -5.18 . (78) 



Mpi 

The fixed lattice case, Eq. (77b) is consistent with Eq. ([78]) only if N c \ < 8.96, which would suggest that loop quantum 
gravity effects should be observable in the spectrum of the CMB temperature anisotropies. In the lattice refinement 
case however no such incompatibility between theoretical and observational constraint exists. 

In conclusion, lattice refinement is required to naturally obtain a successful inflationary model. In contrast to 
this the fixed lattice case requires the number of efoldings of inflation that can be considered classical to be severely 
fine-tuned to match the observational data, leading to the conclusion that most of the inflationary era requires loop 
quantum gravity corrections to be included. This is at odds with the startling agreement between classical CMB 
calculations and observation. 

Finally, it is important to note that whilst ft = /io/i -1 ^ 2 as been suggested as a natural choice for the scaling 
behaviour of the lattice Q, there is no rigorous justification for this choice. A different choice would result in a 
different Hamiltonian constraint equation, Eq. (|52p , and all the subsequent analysis would change. Dynamic lattice 
refinement represents a large freedom in the formulation of LQC that, at present, is fixed only by phenomenological 
arguments. However, deriving the form of this refinement from the exact LQG equations remains, at present, not 
possible. 

At last, but not least, we would like to bring to the attention of the reader that an alternative approach to the one 
we followed here, is to perform a stability analysis of the dynamical difference equations, as suggested in Ref. [9(. More 
precisely, to investigating the scale at which the lattice is unable to support the continuum oscillations, one could 
perform a von Neumann analysis of the difference equations on the grid [2(| ■ Unfortunately, applying this method to 
a refining lattice case is more complicated since plane waves could no longer be used. 
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